LATTICE CALCULATION OF THE GLUON PROPAGATOR AND 
THE STRANGENESS MAGNETIC MOMENT OF THE NUCLEON 



Anthony G. Williams 
CSSM and Department of Physics and Mathematical Physics 
University of Adelaide, Australia 5005 
awilliamOphysics . adelaide . edu.au 
tittp : //www . physics . adelaide . edu . au/ cssm 



ABSTRACT 

This contribution summarizes a presentation combining two topics in lattice 
QCD. Firstly, the gluon propagator in Landau gauge is calculated in quenched 
QCD on a large (32 3 x 64) lattice at (3 = 6.0. New structure seen in the infrared 
region survives conservative cuts to the lattice data, and serves to exclude a 
number of models that have appeared in the literature. Secondly, I report on 
a recent lattice QCD calculation of the strangeness magnetic moment of the 
nucleon. The result is Gj^(0) = —0.36 ± 0.20. The strangeness Sachs electric 
mean-square radius {rfjE was also found to be small and negative. 



1. Introduction 

Here I summarize the central results described in a single long parallel session 
talk, which was comprised of two distinct studies of lattice QCD. The first of these 
is concerned with the direct lattice calculation of the Landua gauge nonperturbative 
gluon propagator, an understanding of which is central to our understanding of the 
nature of confinement. The second topic is a lattice calculation of the strangeness 
magnetic moment of the nucleon. Different confining quark models predict a variety of 
results for this quantity and so lattice calculations provide a benchmark against which 
we can test our understanding of hadron structure. 

2. Lattice Calculation of the Gluon Propagator 

2.1. Motivation 

The infrared behaviour of the gluon propagator is important for an understanding 
of confinement. Previous coiuBctures range from a strong divergenceafi to a propagator 
that vanishes in the infrared.™ Lattice QCD should in principle be able to resolve this 
issue by first-principles, modebindependent calculations. However, lattice studies have 
been inconclusive up to now ,00 since they have not been able to access sufficiently 
low momenta. The lower limit of the available momenta on the lattice is given by 



Oram — 27r/L, where L is the length of the lattice. Here we will report results using a 
lattice with a length of 3.3 fm in the spatial directions and 6.7 fm in the time direction. 
This gives us access to momenta as small as 400 MeV. 

2.2. Summary of Formalism 

The gluon field can be extracted from the link variables U^(x) using U^x) = 
e igoaA^(x+jx/2) _|_ q(q%^ Inverting and Fourier transforming this, we obtain 



\{UM ~ Ul(-q)) ~ gTr (UM - UH-q) 



(1) 



where U^(q) = J2x e ~^ x U^(x) and A^(q) = t a A^(q). The available momentum val- 
ues q are given by q^ = 27rn^/ (aN^), = 0, ... , — 1, where A^ is the number 
of points in the /i direction. The gluon propagator D^ u (q) is defined as D® b u (q) = 
(A^(q)A b u (—q)) /V. The Landau gauge propagator in the continuum has the struc- 
ture D ab u (q) = 5 ab [Sfj, v — {q^qv/q 2 )]D(q 2 ). At tree level, D(q 2 ) will have the form 
D^°\q 2 ) = 1/q 2 . On the lattice, this becomes 



£ (0) (<7) = 1/E -sin^H • (2) 




Since QCD is asymptotically free, we expect that up to logarithmic corrections, q 2 D(q 2 ) — > 
1 in the ultraviolet. Hence we define the new momentum variable q by q^ = (2/a) sin(g M a/2) 
and work with this throughout. The (bare) lattice gluon propagator is related to the 
renormalised continuum propagator D R (q;fi) via D L (qa) = Z^(fi, a)D R (q; fi). The 
renormalisation constant Z 3 (fi,a) can be found by imposing a momentum subtrac- 
tion renormalisation condition Df>(q)\ q 2 =fl 2 = 1/fi 2 . The asymptotic behaviour of the 
renormalised gluon propagator in the continuum is given to one-loop level by D R (q) = 

D haie (q)/Z 3 = (1/q 2 ) [(1/2) \n(q 2 / A 2 )}~ dD with d D = (39 - f - 4JV>)/[4(33 - 2N f )} = 
13/44, where both the gauge parameter £ and the number of fermion flavours Nf are 
zero in this calculation. 

2.3. Simulation Parameters, Finite Size Effects and Anisotropies 

We have analysed three lattices, with different values for the volume and lattice 
spacing. The details are given in Table [I]. In the following, we are particularly interested 
in the deviation of the gluon propagator from the tree level form. We will therefore 
factor out the tree level behaviour and plot q 2 D(q 2 ) rather than D(q 2 ) itself. 

Fig. [I] shows the gluon propagator as a function of qa for the small and large 
lattices, with momenta in different directions plotted separately. For low momentum 
values on the small lattice, there are large discrepancies due to finite size effects between 
points representing momenta along the time axis and those representing momenta along 
the spatial axes. These discrepancies are absent from the data from the large lattice, 
indicating that finite size effects here are under control. 



Name (3 a st (GeV) Volume iV conf gg M a X 



Small 


6.0 


1.885 


16 3 x 48 


125 


2tt/4 


< io~ 6 


Large 


6.0 


1.885 


32 3 x 64 


75 


2tt/4 


< 10~ 6 


Fine 


6.2 


2.63 


24 3 x 48 


223 


2tt/4 


< 10" 6 



Table 1. Simulation parameters 
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Fig. 1. Componentwise data, for the small lattice (left) and the large lattice (right). The 
filled triangles denote momenta directed along the time axis, while the filled squares denote 
momenta directed along one of the spatial axes. 



However, at higher momenta, there are anisotropics which remain for the large 
lattice data, and which are of approximately the same magnitude for the two lattices. In 
order to eliminate these anisotropies, which arise from large momenta (i.e., finite lattice 
spacing errors), we select momenta lying within a cylinder of radius Aqa = 2 x 27r/32 
along the 4-dimensional diagonals in order to distribute large momenta across the four 
momentum components.L 

2.4- Scaling behaviour 

Since the renormalised propagator D R (q; /z) is independent of the lattice spacing, 
we can derive a simple, g-mdependent expression for the ratio of the unrenormalised 
lattice gluon propagators at the same value of q: 



D c (qa c ) Z 3 (fi,a c )D R (q; fi)/a 2 c 



D f (qa f ) Z 3 (/j } a f )D R (q;n)/a 



f 



Z f a? c 



(3) 



where the subscript / denotes the finer lattice (/3 = 6.2 in this study) and the subscript 
c denotes the coarser lattice (@ = 6.0). We can use this relation to study directly 
the scaling properties of the lattice gluon propagator by matching the data for the 
two values of (3. This matching can be performed by adjusting the values for the 
ratios Rz = Zj/Z c and R a = aj/a c until the two sets of data lie on the same curve. 
Fig. shows the data for both lattice spacings as a function of qa before shifting. This 




Fig. 2. The dimensionless, unrenormalised gluon propagator as a function of ln(ga) for the 
two values of (5. The triangles denote the data for the small lattice at /3 = 6.0, while the 
circles denote the data for (5 = 6.2. 



procedure gives us an estimate of R a = 0.745i| 2 , which is fully compatible with the 
string-tension valuel of 0.716 ± 0.040. The corresponding estimate for the ratio of the 
renormalisation constants is Rz = 1.038i|i- 

2.5. Model Fits 

We have demonstrated scaling in our lattice data over the entire range of q 2 
considered, and will now proceed with model fits. The following functional forms have 
been considered: 



D(q 2 ) = 

D(q 2 ) = 
D(q 2 ) = 

D(q 2 ) = 

where 

D(q 2 ) = 
D(q 2 ) = 
D(q 2 ) = 
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(5) 
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(8) 
(9) 
(10) 



We have also considered^ special cases of the three forms d8|)- (pTJ|) , with M IR = 
M uv or with specific values for the exponent a. Equations (f|) and ([5]) are modified in 
order to exhibit the appropriate one-loop asymptotic behaviour. Models (§) and @ 
are constructed as generalisations of (Jm with the correct dimension and asymptotic 
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0.23t\ 4 






1 


1998 


2.2 





0.23 






1 1 


163 


2.41 + x ° 2 




0.14+A 




0.2911 


l_J 


50.3 


6.511 




0.24 +x 3 6 


0.271? 






2.96 


1.54+g 


i.24±ii 


0.46+u 


0.96lf 7 


1.311,1 


g M IR = M UV 


3.73 


1-71±8 


0.84+ 2 °9 


0.48 +i 2 7 




1.521P 


I 


1.57 


1.78l| 


0.49t6 7 


0.43lf 


0.201? 7 , 


0.9511 


g M IR = M UV 


4.00 


1.6211 


0.58±f 


0.40 1^ 




0.92t\ 7 




47 


2.091?!] 


29+ 166 


0.22+i° 6 


0.14±iS 


0.49_ + i° 6 



Table 2. Parameter values for fits to models (|i|)-(|To|). The values quoted are for fits to the en- 
tire set of data. The errors denote the uncertainties in the last digit (s) of the parameter values 
which result from varying the fitting range. The fitting ranges considered when evaluating 
the uncertainties are those with a minimum of 40 points included and with the minimum 
value for qa no larger than 1.03 (point number 40). 



behaviour. All models were fitted to the large lattice data using the cylindrical cut. 
The lowest momentum value was excluded, as the volume dependence of this point 
could not be assessed. In order to balance the sensitivity of the fit between the higli- 
and low-momentum region, nearby data points within A(qa) < 0.05 were averaged. \ 2 
per degree of freedom and parameter values for fits to all these models are shown in 
Table |2|. It is clear that model (|^) accounts for the data better than any of the other 
models. The best fit to this model is illustrated in Fig. |3[ 

2. 6. Summary of Numerical Results 

We have evaluated the gluon propagator on an asymmetric lattice with a large 
physical volume. By studying the anisotropies in the data, and comparing the data 
with those from a smaller lattice, we have been able to conclude that finite size effects 
are under control on the large lattice. A clear turnover in the behaviour of q 2 D(q 2 ) 
has been observed at q ~ lGeV, indicating that the gluon propagator diverges less 
rapidly than 1/q 2 in the infrared, and may be infrared finite or vanishing. The data 
are consistent with a functional form D(q 2 ) = Di^ + -Duv, where 

M ~ 1 GeV, and -Duv is the appropriate asymptotic form. A more detailed analysis^ of 
the asymptotic behaviour reveals that the one-loop formula remains insufficient at q 2 = 
50GeV 2 . Issues for future study include the effect of Gribov copies and of dynamical 
fermions. We also hope to use improved actions to perform realistic simulations at 
larger lattice spacings. This would enable us to evaluate the gluon propagator on larger 
physical volumes, giving access to lower momentum values. 




Fig. 3. The gluon propagator multiplied by q 2 , with nearby points averaged. The line illus- 
trates our best fit to the form defined in @. The fit is performed over all points shown, 
excluding the one at the lowest momentum value, which may be sensitive to the finite volume 
of the lattice. The scale is taken from the string tensionJa 

3. Nucleon Strangeness Magnetic Moment 

3.1. Motivation 

We summarize some recent results0 on nucleon electromagentic form factors, in- 
cluding the strangeness electric and magnetic form factors. The strangeness content of 
the nucleon has been a topic of considerable recent interest for a variety of reasons. 
The studies of nucleon spin structure functions in polarized deep inelastic scattering 
experiments at CERN and SLACJ13 combined with neutron and hyperon (3 decays, have 
turned up a surprisingly large and negative polarization from the strange quark. In ad- 
dition, there is a well-known long-standing discrepancy between the pion-nucleon sigma 
term extracted irom the low energy pion-nucleon scattering^ and that from the octect 
baryon masses.E-3 This discrepancy can be reconciled if a significant ss content in the 
nucleon0E3 is admitted. To address some of these issues, an experiment to measure the 
neutral weak magnetic form factor Gfj via elastic parity-violating electron scattering at 
backward angles was recently carried out by the SAMPLE collaboration.^ The strange- 
ness magnetic form factor is obtained by subtracting out the nucleon magnetic form fac- 
tors G P M and G n M . The reported value is G S M (Q 2 = O.lGeV 2 ) = +0.23±0.37±0.15±0.19 
and does not yet provide a strong constraint on the sign. 

Theoretical predictions of G S M (0) vary widely. The values from various models 
and analyses range from— 0.75 ± 0.30 in a QCD equalities analysis^ to +0.37 m an 
SU (3) chiral bag model.ll3 While a few give positive values,t§EJj most model predictions 
are negative with a ty_oical range of —0.25 to —0.45. Summaries of these predictions 
can be found in Refs.E-3EJ A similar situation exists for the strangeness electric mean- 
square radius (t 2 )e- A number of the predictions are positive while a few are negative^ 
Elastic ep and e A He parity- violation experiments are currently planned at TJNAFEJ 
to measure the asymmetry Alr at forward angles to extract (t 2 )e- Hopefully, they will 
settle the issue of its sign. 



3.2. Numerical Results 



The lattice formulatiim. of the electromagnetic and other form factors has been 
given in detail in the past.E3'E3 Here, we shall concentrate on the DI contribution, where 
the strangeness current contributes. In the Euclidean formulation, the Sachs EM form 
factors can be obtained by the combination of two- and three-point functions 

G% a N (t,p) =T / e-^(0\ X a (x)x a (0)\0) (12) 

X 

G a N ^ N (t f ,p,t,q) = £ e-^ + ^(0| X a (o; / )^(x)x /3 (0)|0), (13) 

Xf,X 

where \ a is the nucleon interpolating field and V^(x) the vector current. With large 
Euclidean time separation, i.e. tf — t » a and t » a, where a is the lattice spacing, 

TtG a N % N (tf, 0, t, q) G%« N (t, 0) e ijk q k „ . 2 



Gff N {t f ,0) G$? N (t,q) E q + m 

T^G%y^^{tf, 0, t, q) G% a N (t, 0) ^ 

G% a N (t f ,0) G% a N (t,q) 



Gm(q ), (14) 
G E (q 2 ), (15) 



where Tj = cr.j(l + r y^)/2 and = (1 + 7 4 )/2. 

We shall use the conserved current from the Wilson action which, being point- 
split, yields slight variations on the above forms and these are given in Ref.£3 Our 50 
quenched gauge configurations were generated on a 16 3 x 24 lattice at (3 = 6.0. In the 
time direction, fixed boundary conditions were imposed on the quarks to provide larger 
time separations than available with periodic boundary conditions. We also averaged 
over the directions of equivalent lattice momenta in each configuration; this has the 
desirable effect of reducing error bars. Numerical details of this procedure are given in 
Refs.BBJ The dimensionless nucleon masses M^a for k = 0.154, 0.152, and 0.148 are 
0.738(16), 0.882(12), and 1.15(1) respectively. The corresponding dimensionless pion 
masses m^a are 0.376(6), 0.486(5), and 0.679(4). Extrapolating the nucleon and pion 
masses to the chiral limit we determine k c = 0.1567(1) and rriMa = 0.547(14). Using 
the nucleon mass to set the scale to study nucleon properties,c3'cJ the lattice spacing 
a -1 = 1.72(4) GeV is determined. The three k's then correspond to quark masses of 
about 120, 200, and 360 MeV respectively. 

The strangeness current s7 M s contribution appears in the DI only. The full details 
of the extraction can be found in Ref.EJ and we satisfy ouselves here by simply quoting 
the results. We use a monopole form to extrapolate G s M (q 2 ) with nonzero q 2 to G S M (0), 
giving G S M (0) = —0.36 ±0.20. Correlations are taken into account and the error is from 
jackknifing the fitting procedure. A similar analysis is done for the strange Sachs electric 
form factor G s E (q 2 ) and we find that G S E (0) is consistent with zero as it should be. We 

find for the electric mean-square radius {t 2 s )e = —6—^T^-\q 2 =o — —0.061 ± 0.003 fm 2 . 
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